Electronic properties of the armchair graphene nanoribbon. 
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We investigate the electronic band structure of an undoped graphene armchair nanoribbon. We 
demonstrate that such nanoribbon always has a gap in its electronic spectrum. Indeed, even in 
the situations where simple single-electron calculations predict a metallic dispersion, the system is 
unstable with respect to the deformation of the carbon-carbon bonds dangling at the edges of the 
armchair nanoribbon. The edge bonds' deformation couples electron and hole states with equal 
momentum. This coupling opens a gap at the Fermi level. In a realistic sample, however, it is 
unlikely that this instability could be observed in its pure form. Namely, since chemical properties 
of the dangling carbon atoms are different from chemical properties of the atoms inside the sample 
(for example, the atoms at the edge have only two neighbours, besides additional non-carbon atoms 
might be attached to passivate unpaired covalent carbon bonds), it is very probable that the bonds 
at the edge are deformed due to chemical interactions. This chemically-induced modification of the 
nanoribbon's edges can be viewed as an effective field biasing our predicted instability in a particular 
direction. Yet by disordering this field (e.g., through random substitution of the radicals attached 
to the edges) we may tune the system back to the critical regime and vary the electronic properties 
of the system. For example, we show that electrical transport through a nanoribbon is strongly 
affected by such disorder. 
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I. INTRODUCTION 

Graphene is attracting considerable attention due to 
its unusual electronic properties including: large mean 
free path, "relativistic" dispersion of the low-lying elec- 
tron states, and "valley" degeneracy [1]. These remark- 
able features made many researchers hope that some 
day graphene mesoscopic structures might revolutionise 
nanoscience. Thus, a substantial amount of effort has 
been invested investigating graphene devices, such as 
quantum dots [2] , bilayer structures [3] , and nanoribbons 

[4]- 

Studying the physics of nanoribbons, a certain discrep- 
ancy between results of first-principle calculations [5, 6], 
experiments [7, 8], and single-electron approximations 
[9, 10, 11, 12] was stumbled upon: whereas the single- 
electron approximation predicts that, depending on its 
width, an armchair nanoribbon could be either semicon- 
ducting or metallic, both experiments and first-principle 
calculations suggest that it is always semiconducting: a 
nanoribbon of any width W has gap at the Fermi level; 
the magnitude of the gap scales as 1/W. 

Two different mechanisms were proposed to explain 
this disagreement. According to the first mechanism the 
gap is due to electron-electron interactions [13]. At zero 
doping the interaction induces a charge gap. The physics 
here is similar to the Mott transition in the Hubbard 
model at half filling. If the Coulomb coupling constant 
g is large (g > t, where t is the carbon-carbon hopping 



amplitude), the gap scales as g. If g is small, the gap 
vanishes faster than g. Since g cx 1/W, such mechanism 
is consistent with the observed scaling A cx 1/ W for small 
W only. For large W the gap decays faster than 1/W. 

Another way to explain the gap was outlined in [5], 
where first-principle computations have shown that the 
length of the carbon-carbon bonds at the hydrogen- 
passivated edge is shorter than the length of the bonds in 
the bulk. Due to this, the hopping amplitude across the 
bonds dangling from the nanoribbon edge (i e dge) differs 
from t: t c d g0 = t + 6t, St ~ O.lt (see Fig.l where the 
short bonds are shown in bold). Because of St, the na- 
noribbon's Hamiltonian acquires an additional term. It 
couples states with the same momentum above and be- 
low the Fermi energy; such coupling opens a gap in the 
electronic spectrum. Assuming that St depends mainly 
on the chemical properties of the edge (e.g., nature of the 
passivating radical or absence thereof), one can demon- 
strate that the effective strength of this coupling is in- 
versely proportional to W. Thus, the relation A oc 1/W 
is recovered. 

In addition to these mechanisms, other instabilities can 
be present as well. Which mechanism of gap generation 
dominates depends on a variety of factors, such as electro- 
static screening by the gate electrode, mechanical forces 
applied to the nanoribbon, edge passivation, etc. 



A. Summary of our results 

In this paper we will study in detail the second mech- 
anism described above. We will show that the armchair 
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nanoribbon is unstable toward the deformations of the 
edge bonds. Specifically, we prove the following: if the 
radical passivating the edge is carefully picked (to en- 
sure that it docs not modify the hopping at the edges: 
Wgc = t) , then the total energy of the nanoribbon would 
become a decreasing function of the edge bond deforma- 
tion. In other words, spontaneous deformations of the 
dangling bonds are favourable for they decrease the total 
nanoribbon energy. Such deformations make i d go un- 
equal to t. 

To understand the physics behind this, it is convenient 
to describe such deformation using a real one-component 
(Ising-like) order parameter St. At the mean-field level 
we can state that the energy of the deformed bonds is 
increased by the amount £b ~ (St) 2 > 0, while the energy 
of the conducting electrons decreases by the amount e e i ~ 
A 2 ln(A/t) < 0. Since A — \5t\, the total energy, eb + e e i, 
always has a minimum at not-zero \St\. 

It is useful to extend the analogy between our system 
and the Ising model further. Notice that our chemical 
forces deforming the edge bonds are similar to the "mag- 
netic field" coupled to the order parameter (the "mag- 
netisation") in the Ising model. Thus, it is natural to 
refer to these chemical forces as the "edge field" . 

In practice, we have very limited control over the mag- 
nitude of the "edge field" : it is very difficult to select the 
radical attached to the nanoribbon's edge to guarantee 
that 5t = 0. Using the analogy between our nanoribbon 
and the Ising model, we conclude that, since the "field" 
is always switched on, our system is always away from 
the critical regime. 

Yet even a system with a strong "edge field" may be 
tuned to criticality by disordering the "field" . This way, 
the disorder may be used as a tool to control spectral 
properties of nanoribbons. 

Disorder can be introduced by passivating the edges 
by radicals of two different types, randomly distributed 
along the length of the nanoribbon. The hopping ampli- 
tude at the edge i e dge would become a function of the 
coordinate along the ribbon: 

W g c (x) = t + 5t + 5t dis (x), ( 1 ) 

where St is independent of x, while Stdi s (x) is the disor- 
dered part. 

In such system the electronic spectrum is very sensitive 
to the relative strengths of St and St^is'- when the disorder 
is weak, the spectrum has (pseudo)gap due to non-zero 
St term; when St = 0, the pscudogap closes, and the 
nanoribbon is in the disorder-dominated regime. 

The crossover between these two regimes can be ob- 
served with the help of transport measurements. In a 
nanoribbon whose length L is thermodynamically large, 
the conductance vanishes exponentially in both regimes: 
when there is no disorder, the finite gap removes density 
of states from the Fermi level; when there is disorder, 
the wave functions are localized. Thus, for very long 
nanoribbons the conductance is suppressed regardless of 
the disorder strength. 



The situation is different for a mcsoscopic sample. 
Note that the localization is a weaker phenomenon than 
the gap opening in the sense that the inverse of the lo- 
calization length is the second order in the disor- 
der strength: = 0(<5t^ is ), while the inverse of the 
length scale /~ a p ~ A/wp charachterizing the gap size is 
linear in St: — 0(\St\). That is, a longer nanorib- 
bon is required to observe a well-developed localization. 
Therefore, it is possible to choose L sufficiently short to 
have no localization features at any disorder, and yet 
long enough to observe the spectral gap when the disor- 
der is low. Thus, when the disorder is absent, the low- 
temperature conductance of the sample is exponentially 
suppressed due to the gap; otherwise, the conductance is 
finite down to the lowest temperatures. Therefore, the 
disorder effectively closes the gap. 

The paper is organized as follows. In Sect. II we derive 
the model for a nanoribbon with deformed edges. The 
instability of the nanoribbon toward the deformation of 
the edge bonds is discussed in Sect. III. In Sect. IV we 
investigate the effect of disorder on the transport proper- 
ties of the nanoribbon. Sect. V presents the conclusions. 



II. MODEL 

In this section we will obtain the Hamiltonian for a 
graphene nanoribbon with deformed edge bonds. Our 
derivation relies on basic facts of the graphene physics 
which are discussed in Ref. [1]. 

A. Tight-binding model for graphene 

For completeness, in this subsection we quickly red- 
erive basic single-electron properties of a graphene sheet. 
This gives us an opportunity to introduce notation we 
will need below. 

It is common to describe a graphene sample in terms 
of a tight-binding model on a honeycomb lattice. Such 
lattice can be split into two sublattices, denoted by A 
and B. 

The Hamiltonian of a graphene sheet is given by: 

H = -tJ2 E c rCr +5i +H.c, (2) 

R i=l,2,3 

where R runs over sublatticc A. The vectors Si (i — 

1,2,3) connect the nearest neighbours. These are (see 
Fig. 3): 

5 1 = oo(-l,0), (3) 

5 2 = a (l/2,V3/2) (4) 

5 3 = o (l/2, -V3/2). (5) 

The symbol a denotes the carbon-carbon bond length, 
which is about 1.4 A. 
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The corresponding Schrodinger equation can be writ- 
ten as: 



R+Si+a^ 



i=l,2 



^+S 1 = -t^R-t J2 ^R-«M> 
i=l,2 



(6) 
(7) 



where -0^ (^>- B 



R+<5i 



denotes the wave function value at 



the site R (at the site R + Si) of sublattice A (sublattice 
B). The primitive vectors of the honeycomb lattice are: 



ai = oo(3/2,\/3/2), 
a 2 = oo(3/2, -V3/2). 



(8) 
(9) 



They connect nearest neighbours on the same sublattice. 

Since the primitive cell contains two atoms, it is conve- 
nient to define a two-component (spinor) wave function: 



0: 



R+<5i 



(10) 



With this notation the action of H on spinor \?k can be 
expressed as: 



tk = t 



1 + 2exp 



-tk 

-*£ o 

'V3 



cos 



k y a 



For every k there are two eigenstates: 



1 



exp (z6» k ) 
with the eigenvalues: 



(11) 
(12) 



(13) 
(14) 



e k ± = ±|ik| = ±i 



1 + 4 cos 



3 i. 



cos 



— k y a 



4 cos 2 



(15) 



The states with negative (positive) energy are filled 
(empty) at T = 0. 

The quantity £k vanishes at six points within the 
Brillouin zone (see Fig. 2): (0, ±47r/(3v / 3ao)) and 
(±27r/(3ao),±2vr/(3^a )). These are the locations of 
the famous Dirac cones of graphene. 

These six cones can be split into two equivalence 
classes: the locations of any two cones inside the same 
equivalence class differ by a reciprocal lattice vector. 
Thus, we do not need all six of them. Two inequivalent 
cones are sufficient: 

cone K :\l k = (27r/(3ao),27r/(3\/3a )); (16) 
cone Kf : = (0, 47r/(3V3a )). (17) 



B. The nanoribbon spectrum 

Consider now the electron states of the armchair na- 
noribbon (see Fig. 3). Such nanoribbon is defined by 
the condition < y < W. The width of the nanoribbon 
W is a multiple of \/3ao/2: 



W 



where M is an integer. 



M, 



(18) 



The Schrodinger equation for the edge sites differs from 
Eqs. (6) and (7). For sites on the upper edge (y — W) 
we can write: 



^R+S, 



(19) 
(20) 



Note the absence of the summation over the lattice vec- 
tors in the right-hand side of these equations. This 
boundary condition is called free. It can be easily gener- 
alized for the lower edge. 

An armchair nanoribbon is invariant under a shift over 
3a along the x-axis. Thus, a 3a long nanoribbon seg- 
ment can be thought of as a nanoribbon unit cell. There 
are (M + 1) graphene unit cells in a nanoribbon unit cell. 
The nanoribbon Brillouin zone is: 



77 k 7F 
3ao x 3<2o 



(21) 



Now we assume that all carbon-carbon bonds of our 
system have the same hopping amplitude t (that is, the 
deformation of the edge bonds is absent) . The wave func- 
tion vPr, satisfies the free boundary conditions Eqs. (19) 
and (20) at R = (x, W) and similar conditions at R = 
(x,0). It is not convenient, however, to work with such 
boundary conditions directly. Fortunately, if we add an 
additional row of lattice sites at each edge and demand 
that the wave function vanishes at these auxiliary sites 
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(see Fig. 3), then the wave function at the real edge sites 
satisfies the free boundary condition. In other words, 
our free boundary condition problem for a nanoribbon 
of width W is equivalent to the zero boundary condition 
problem for a nanoribbon of width W + \/3ao- Thus, we 
want: 



*R|y=-V3a /2 



R 



y=W+V3a /2 



= 0. 



(22) 



The eigenfunction of the Hamiltonian H, Eq. (2), satis- 
fying Eq. (22) can be written as follows: 



"(23) 



where c\_i are complex coefficients. Note that this eigen- 
function has a well-defined value of the momentum k x 
along the x-axis, but not of the momentum k y along the 
y-axis, since our system has no translational invariance 
in the y-direction. 

The values of k y and Ci j2 in Eq. (23) must be chosen to 
satisfy Eq. (22). Since the spinor \l/k± remains the same 
when the sign of k y changes, we derive: 



sin 



k y (W + V3, 



Cl 



-exp 



(^-iVSkyao^ 



0, 



C2- 



(24) 
(25) 



We then obtain the following quantization condition: 

2wn 



k = 
V V3(M + 2)a ' 



(26) 



where n is an integer. Thus, the nanoribbon spectrum 
£n(k x ) consists of a set of one-dimensional branches la- 
belled by an integer, n. 

Using the above results it is possible to construct ex- 
plicitly the nanoribbon eigenfunction. However, for our 
purposes it is more convenient to define an effective 
Hamiltonian for a given branch. Let us look for a na- 
noribbon Hamiltonian eigenfunction in the form: 



*r = ^n{x) sin 



27m 



V3(M + 2)a \ 
where the spinor ^f(x) is defined as: 



y + 



V3ao 



, (27) 



*n(*) = ,/JJ 



ip%(x-a ) 



(28) 



In this equation x = 3a m/2, and m is an integer. 

Substituting the expression for \I/ n into the Schrodinger 
equations (6) and (7), we obtain: 





-t 



V n (x) + 2cos(V3k y a /2) 



*„(x + 3a /2) + 




-t 



*„(a; - 3a /2) 



(29) 



In fc-space this equation has the form e^ n = H(n)$f n , 
where the effective Hamiltonian for branch n is: 



H(n) 





-/* 







1 + 2 cos 



whi) cxp (-i kxao 



, (30) 
• (31) 



This Hamiltonian possesses accidental symmetry: the na- 
noribbon remains unchanged under a shift by 3ao along 
the x-axis, yet the effective Hamiltonian H(n) is invariant 
under a shift by 3ao/2. That is, the effective symmetry of 
the Hamiltonian is higher than the geometric symmetry 
of the underlying system. This symmetry is destroyed 
when the edge bonds are deformed. However, since in 
this paper we are interested in the low-energy properties 
of the system, this peculiarity will play no role in what 
follows. 

Within our labelling scheme the same branch may ap- 
pear under different values of index n. Obviously, n 
and — n correspond to the same branch. Furthermore, 
if n > 0, n' > 0, and n — — n'mod (M + 2), both n and 
nf define the same branch. This means that there are 



(M + 1) independent branches: 

< n< M + 2. 



(32) 



This is precisely the number of graphene unit cells in a 
nanoribbon unit cell. The Hamiltonian of the nanoribbon 
is the direct sum of the H(n)'s: H = J2 n H(n). 

The dispersion associated with a specific branch is 
given by Eq. (15), where k y is fixed by Eq. (26). 
In other words, a branch samples the function along 
the line k y = const. 

Let us now find out under which circumstances the 
gapless branches appear. Branch n [n satisfies Eq. (32)] 
is gapless if the complex equation (a system of two real 
equations) tfc x ,„ = has a root satisfying Eq.(21). 

Solving this system of trigonometric equations one can 
prove that a zero eigenvalue at k x = appears when 
cos(t™/(M + 2)) = -1/2. The latter condition is fulfilled 
if the argument of the cosine is 27r/3. This is possible 
provided that n — no, where: 



n 



2(M + 2) 



(33) 



and no is an integer. In other words, the gapless branch 
is present only when (M + 2) is divisible by 3. This 
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agrees with the conclusions of Rcf. [12], where the single- 
electron calculations for a nanoribbon with no edge de- 
formation was performed. 

The effective Hamiltonian for the gapless branch is: 



we find: 



Ho = H(no) 







1 — exp — i 





.3k x a 



e k=c± = ±2t 



■,111 | -k x a 



(34) 
(35) 
(36) 



The gapless branch is characterized by k y = 47r/(3\/3ao)- 
Since k = (0, 47r/(3\/3ao)) is the location of the cone JC' 
[see Eq. (17)], one can say that the gapless branch is 
found in the nanoribbon spectrum only when the quan- 
tization condition Eq. (26) allows for existence of the 
branch passing through the cone JC' . 



III. SPONTANEOUS GENERATION OF THE 
GAP 

In this section we show that the electronic branch, 
which appears to be gapless according to the calculations 
reported above, is, in fact, unstable toward the sponta- 
neous opening of the gap. We prove that the edge bond 
deformation is one possible instability leading to the gap 
generation. 



Modification of the Hamiltonian due to edge 
deformation 



To establish such an instability we need to calculate the 
ground state energy of the nanoribbon with edge bonds 
deformed as shown on Fig. 1. 

To achieve this aim we first determine how the edge de- 
formation affects the Hamiltonian of the gapless branch. 
We denote by SH = SH\ + SH^ the contribution to the 
Hamiltonian due to edge deformation. The subscript 
T ('h') corresponds to a bond deformation at the lower 
(higher) edge of the nanoribbon (see Fig. 3). 

The matrix element of 8H\ between two states is equal 
to: 



($\6H 1 \M>) = -5tJ2KK m +c-C- 



(37) 



The summation in this formula runs over the deformed 
bonds at the lower edge: R m = (3aom, 0). 

Substituting in Eq. (37) wave functions consistent with 
Eq. (22), i.e., 



R 



*e 



sin 



27m 



V3(M + 2)a 



IJ + 



V3< 



ao 



(38) 



<Je 



z h.- j- .x 



2tt s 



sm 



( 47ry 

VSv^ao 3 / 



|*> = * f ( ° Xf nM *sin 2 (2V3) (39) 



-St 
x exp [— 3iao{k x — k' x )m] . 

m 

Therefore, for k x ,k' x satisfying Eq.(21), Hamiltonian 8H\ 
equals to: 



SHi 



3<5i / 1 



8M + 8 V 1 °/ K ' K 



(40) 



where 2(M + 1) in the denominator comes from the wave 
function normalization. 

It is trivial to demonstrate that SH\ = SH^. Superfi- 
cially, this identity appears to be incorrect: clearly, there 
should be a difference between the lower and higher edges 
at the level of the Hamiltonian. However, we must re- 
member that the expression Eq. (40) is not the total 
Hamiltonian 8H\ (which is indeed different from SH^), 
but rather its projection on the subspace spanned by a 
specific branch. These projections cannot discriminate 
between the lower and the higher edge. 

The total edge Hamiltonian SH is equal to twice 5H\. 
Thus, the Hamiltonian for a nanoribbon with the de- 
formed edge bonds is: 



H + 5H 





-5t e f[ — tt 



St e S — tk x 



3St 



8t e S = 



(41) 
(42) 



4M + 4' 

whose dispersion is given by £fc x ± = ±e(k x ), where: 

s(k x ) = {5t cS f + 4t(t + St e ff)sin (3fca;ao/4),(43) 

which has a gap A = 2 1 <5t e fT ] - 

Since our formerly gapless branch now acquired the 
gap, it might be confusing to refer to such branch as 
'gapless'. Instead, we will now call it 'no branch', where 
n is given by Eq. (33). 

Heuristically, one can say that the gap appears because 
the boundary conditions at the edges have changed. In- 
deed, we explained that the gapless branch exists because 
the quantization condition Eq. (26) makes this branch 
pass through the Dirac cone. When the edge structure 
is altered, the boundary conditions are altered as a re- 
sult. The latter induces a modification of the quantiza- 
tion rule. Thus, in fc-space, the no branch shifts slightly 
off the cone pinnacle's location and acquires a gap. 



B. Edge instability 

To demonstrate the existence of an edge- induced insta- 
bility we need to calculate the energy of the no branch: 



t/(3io) AU 

<k x )^. 

-tt/(3o ) 



2tt 



(44) 
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It is easy to show that: 

2(^eff) 2 



e el /L « e°JL - 



In 



Steff 



+ const. 



(45) 



where e^j is the ground state energy calculated at St e g = 
0. The next term is the most singular correction to 
due to 5t c g. This correction is not analytic in St c g. 

One also needs an expression for the nanoribbon lattice 
energy due to the bond deformation: 



1 ku 1 

eb/L = 2x x = ku 

3a 2 3a 



(46) 



where u is the variation in the edge bond length, k is the 
stiffness of the bond (the energy eh /L is composed of the 
energy of two deformed bonds per unit cell of the nano- 
ribbon; the energy of a single deformed bond is ku 2 /2). 

To proceed further we need to know how to relate the 
deformation of the bond u and St. Such information 
may be extracted from the quantum chemical calcula- 
tions [14]. However, to demonstrate that our system is 
unstable it is enough to assume that at small u we have 
St <~ u. Then one can write the following expression for 



e h /L 



AM + lf 2 
C ; (Ot c g) . 



Oq 



(47) 



where £ is a phenomcnological constant. 

Finally, it is straightforward to check that the total 
nanoribbon energy (e^ + e cl + e b ) has two minima at 
8t c g = ±St* , where: 



St* ~ iexp [-3vrfC(M + l) 2 ] ^ 0. 



(48) 



This expression shows that the nanoribbon energy is 
smallest when the bonds at the edges are deformed and 
the no branch has a gap. 

Note that the calculations presented above rely on the 
mean-field approximation. The latter is applicable to 
our one-dimensional system since the order parameter 
(the bond deformation) is Ising-likc. Thus, no Goldstone 
mode is present, and we do not have to worry about 
critical fluctuations. 

Above we demonstrated that an armchair nanoribbon 
of any width has a gap in its electronic spectrum. There 
are two caveats to our discussion, however. 

First, we proved that for our system at least one gap- 
opening instability exists. We did not prove, yet, that 
the discussed mechanism is the only possible path to 
generate the spectral gap. For example, the electron- 
electron interaction can induce a gap [13]. Ultimately, 
the strongest instability must be determined by compar- 
ing the energies associated with particular mechanisms. 
The energy St* characterizes the strength of the edge de- 
formation instability. This energy scale quickly vanishes 
for wider nanoribbons (M is large) or "stiff" edge bonds 
(C is large). Under these conditions other mechanisms 
might be important. 





Ising model 


Nanoribbon 


Order 


Magnetisation 


Bond deformation u, 


parameter 


S 


u ~ 5t 


External field 


Magnetic field H 


"Edge field" / 


Energy 


aS 2 + bS 4 - HS 


au 2 In u + bu 2 — fu 



TABLE I: Analogy between the Ising model and the nano- 
ribbon of graphene. 



Second, due to special properties of the carbon atoms 
at the edge, it is likely that the edge bonds are deformed 
by chemical forces, which are more powerful than any 
intrinsic instability, including the one we have discussed. 
First-principle numerical simulations support this point: 
analysing Fig. 3 of Ref. [5], we note that the edge bonds 
are deformed even in nanoribbons where the unstable no 
branch is absent; moreover, the deformation magnitude 
is independent of the nanoribbon width. 

One can say that our order parameter St cS is coupled to 
the fictitious "field" , which takes the system away from 
criticality (see Table I). In such a situation St c s ^ 0; 
however, the value of St e g and corresponding spectral gap 
is determined not by Eq. (48) but rather by the strength 
of the external "edge field" . 

Thus, the nanoribbon's spectral properties are con- 
trolled by the chemical structure of the edges. If we find 
a way to vary the effects of the chemical edge structure, 
we may tune the electronic properties of the nanoribbon 
to our needs. 

This problem can be dealt with the help of three dif- 
ferent approaches. First, one can try to hand-pick pas- 
sivating radicals to guarantee preservation of the edge 
structure. Second, one can treat different edges of the na- 
noribbon with different radicals, which deform the edge 
bonds in opposite direction. Clearly, these two propos- 
als require considerable experimental work. A third ap- 
proach seems more promising: to close effectively the gap 
it is enough to disorder the "edge field" . We will examine 
this idea in the next section. 



IV. NANORIBBON WITH EDGE DISORDER 

As we have seen in the previous section, edge bond 
deformations can open a gap in the electronic spectrum 
of an armchair nanoribbon. The magnitude of the gap is 
determined by the chemical properties of the passivating 
radical at the nanoribbon's edges. In this section we 
investigate how the chemical disorder at the edges affects 
the nanoribbon's spectrum. We will see that the disorder 
effectively weakens the "edge field" . 

We consider a nanoribbon whose edges are treated by 
two different radicals, 'a' and '/?', which distort the bonds 
in opposite directions: 



St a > 0, 
St/3 < 0, 



(49) 
(50) 
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\St a \ = \St p \ = 5to. 



(51) 



When these radicals randomly attach to the nanoribbon's 
edges, the "edge field" becomes disordered. By adjusting 
the concentrations n a ,/3 of the two radicals, it is possible 
to vary the relative strengths of St and 5t d i s : 

St = n a 5t a + npStp = St (l — 2np), (52) 
St dis {x) = St(x) -Si. (53) 

These two quantities are defined by Eq. (1). 

The effective Hamiltonian for the tiq branch in the 
presence of disorder is equal to: 

H = H + SH{x), (54) 

5H{x) = -^-^ iyTt + St dis (x)}. (55) 

The coordinate representation of 5H(x) may be obtained 
through a procedure similar to the derivation of 5H\ in 
Sect. II. Below we will assume a Gaussian distribution 
law for the random quantity 8t d i s with the correlation 
function: 

(5t dis (x) St dis (x')) = a 2 f((x - x')/a), (56) 
or, in Fourier space: 



Here: 



{St d i s ,k x Sidis-v ) = cr 2 aL f(ak x ) 5k x ,k' x - (57) 



a 2 ~{St f{np-n% (58) 



sets the scale for the disorder strength variation, the scale 
a is the disorder correlation length. The function f(z) is 
a "broadened (5-function" . It is even and non-negative, 
vanishes quickly for \z\ > 1. In addition, this function is 
normalised by the condition: 



/+oo 
-OO 



dzf(z) = 1. 

Its Fourier transform / satisfies: 
/(0) = 1. 



(59) 



(60) 



In the plane-wave basis, the Hamiltonian SH can be writ- 
ten as: 



SH 



1 



AM + 4 V 1 



St 8,, 



*..*; + V"r T *..*i (61) 



=5idi s ,fc x -/s'(62) 



(\T kx , k ,/)=f(a(k x -k' x )l63) 

In Eqs. (61) and (62), the length L appears in the de- 
nominator due to wave function normalization. In these 
equations we defined the dimensionless random field t^m' 



to show explicitly how the disorder matrix elements scale 
with the nanoribbon length L [15]. 

As we see from Eq. (61), the disordered and homoge- 
neous parts of the edge Hamiltonian enter as two differ- 
ent terms. Each term induces a specific modification of 
the nanoribbon's spectrum: the ordered part opens a gap, 
the disorder part localizes the wave functions. The local- 
ization, however, is a weaker phenomenon than the gap 
generation. Intuitively, this sounds quite reasonable: the 
effects of disorder may "average out" to zero, while the 
ordered term acts "coherently" over the whole sample 
length. To make this statement rigorous we will apply 
perturbation theory in orders of SH. We will show that 
for a nanoribbon of a certain length the disorder may 
be treated with the help of perturbation theory, while 
the ordered term may not. To prove this we separately 
consider the two pieces of SH. 



A. Homogeneous edge deformation 5t 

Perturbation theory is applicable when the level spac- 
ing Se ~ ant/L is much larger than the matrix elements 
of SH. 

To establish applicability range of the perturbation 
theory in orders of St, we have to compare Se with the 
gap A. Thus, perturbation theory works if: 



L<&l e&p = a (M + l)(t/6t). 



(64) 



The length scale / gap is analogous to the Compton length 
lc = l/(m e c) of the Dirac equation. 

When the disorder is absent (np = 0), the gap has its 
biggest value, and l gap is the shortest: 



C:* = ao(M + l)(t/St ). 



(65) 



Thus, one can say that perturbation theory in St works 
for any value of np if 



(66) 



B. Disordered edge deformations Stdis 
The perturbation theory in St d - ls is applicable when: 
ant 1 



> 



L M + 1 V L 
This inequality may be transformed into: 

L«a oc = (Af + i) 2 ^4- 



(67) 



(68) 



The scale l\ oc is the localization length. (Our treatment 
of the disordered regime follows closely Ref. [15].) 

When St = (or, equivalently, n a — n@ = 1/2) the dis- 
order is the strongest. We can define the shortest possible 
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localization length ^"" n . It can be estimated as follows. 
If n a — np = 1/2, a has its largest possible value: 



C n = (M + i) 



2 j-2 
an i 



2 "0 



a (6t ) 2 ' 

Therefore, if the sample length satisfies: 
L « C n , 



(69) 
(70) 

(71) 



the perturbation theory in the disorder strength is justi- 
fied for any concentration np. 



C. Conductance of a mesoscopic sample 



It is possible to prove that < ZjJJ™. Indeed, this 



inequality is equivalent to: 



( M + 1 ^ >> f- 



(72) 



Both factors on the left-hand side of this expression are 
much larger than unity. Therefore, unless a is very big, 
'gap ^ 'Joe"- Loosely speaking, this inequality suggests 
that the disordered field is a much weaker perturbation 
than the ordered one. 

Consider now a nanoribbon whose length L satisfies: 



« L « C" 



(73) 



This means that even the strongest disorder (Zi oc = ^"" n ) 
cannot create a well-developed localization in our nano- 
ribbon; on the other hand, when the sample is close to 
perfect order (np <C 1/2) the spectral gap fully mani- 
fests itself. Let us now study the electrical conductance 
of such nanoribbon. 

Assume first that our system has no edge disorder: 
St(x) = St a . Then l gap = Z™, and Eq. (64) is violated. 
Therefore, the "edge field" opens a gap in the spectrum. 
The dimensionless conductance g of a sample with the 
gap is exponentially small at T — 0: 



ln<7 ~ — 



gap 



1 St L 
M + lTcfo' 



(74) 



When we slightly disorder our system by introducing a 
small concentration np of radicals '/?', the conductance 
increases since St decreases [see Eq.(52)]. Thus, the dis- 
order improves electrical conductance! 

In the opposite limit of complete disorder, we have 
St = 0, and the localization length becomes Zi oc = Z™. 
Perturbation theory may now be applied since Eq. (68) 
holds true. Instead of a disorder-induced localization, 
which is a non-pcrturbativc phenomenon, in a sufficiently 
short nanoribbon the disorder creates weak corrections 
to the properties of Hq . In such sample the conductance 
remains finite even at T — 0. 



Thus, we reach the following counter-intuitive conclu- 
sion: a completely ordered nanoribbon shows "insulat- 
ing" behaviour, while a disordered one shows "metallic". 
[We put quotes around "insulating" and "metallic" for 
metal and insulator are quantum phases, which can be 
unambiguously defined only in the thermodynamic limit 
L — > oo; however, the latter limit is incompatible with 
Eq. (73).] 

As the system moves from perfect order to total dis- 
order, it passes through a crossover from "insulating" to 
"metallic" conductance. Indeed, Eq. (74) is applicable 
only when L/l gap is much bigger than unity. If 



L/l 



gap 



1 g ~ 1, 



(75) 



one can validate the perturbation theory not only in or- 
ders of Stdis, but also in orders of St [see Eq.(64)]. There- 
fore, once L exceeds l gap , the exponential dependence of 
g is replaced by a slower function, and g remains of order 
unity down to the completely disordered regime. 



D. Conductance of a long nanoribbon 

Finally, let us comment on the conductance behaviour 
in a thermodynamically large sample, whose length sat- 
isfies: 



'loc 



(76) 



In this situation both terms of 5H cannot always be 
treated with the help of perturbation theory. Conse- 
quently, the conductance is always exponentially sup- 
pressed: when there is perfect order, the conductance 
follows Eq.(74); in the opposite case (complete disorder) 
we have: 



\ng 



(77) 



The latter equation is a manifestation of wave function lo- 
calization, which can only be observed in a sample whose 
length exceeds l\ oc . 

Comparing Eqs. (77) and (74) with the help of 
Eq. (76), we note that, as well as in the case of a short 
nanoribbon, the conductance of a perfectly ordered sam- 
ple is much smaller than the conductance of a totally 
disordered sample. 

The conductance g of a long nanoribbon is a non- 
monotonous function of disorder (see Fig. 4). Such a 
behaviour is a consequence of the crossover from the gap- 
dominated to the disorder-dominated regime. When the 
concentration np is small, the gap is the dominant pa- 
rameter controlling the conductance. Under such cir- 
cumstances Eq.(74) is obeyed. The disorder acts mainly 
to reduce the gap. Thus, if disorder is weak, then g is an 
increasing function of np. 

As np keeps growing, l gap increases, while l\ oc de- 
creases. The crossover occurs at: 



loc 



(78) 
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and the sample enters the disorder-dominated regime. 
The conductance is given by Eq. (77). It is a decreasing 
function of np when the latter is close to 1/2. 

V. CONCLUSIONS 

In this paper we studied the effects of edge bond de- 
formations on the electronic properties of nanoribbons. 
We have seen that a nanoribbon of a certain width is un- 
stable with respect to a spontaneous deformation of the 
edge bonds. While such deformation increases the en- 
ergy of the affected bonds, it also reduces the electronic 
energy. As a result of a 'bond instability' the electronic 
spectrum acquires a gap at the Fermi level. 

We also pointed out that this instability is difficult to 
observe in a real system. The culprit is the chemical 
structure of the edges, which deforms the bonds to opti- 
mize the chemical energy at the edges. Although it might 
be hard to remove this chemical modification of the edge 
bonds, it is quite possible to reduce its effect on the elec- 
tronic spectrum by disordering the radicals passivating 
the edge bonds. We demonstrated that the disorder can 
vary the electrical conductance of a nanoribbon. In case 



of a short nanoribbon, the conductance would change 
from "insulating" regime at low disorder to a "metallic" 
regime at high disorder. When the nanoribbon's length 
is large, the conductance is a non-monotonous function 
of the disorder. Thus, the disorder may be a useful tool 
which allows one to control the electric transport through 
nanoribbons. 
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FIG. 1: A schematic diagram of a graphene sheet, with one 
edge passivated by hydrogen. Deformed bonds are shown by 
bold lines, atoms of carbon and hydrogen are represented by 
the symbols 'C and 'H'. The vectors <5i,2,3 connect nearest 
neighbours on the graphene lattice. The vectors 0,1,2 are 
the primitive lattice vectors. The broken line diamond is the 
graphene unit cell. 
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FIG. 2: Shaded hexagon schematically shows Brillouin zone 
of graphene. The Dirac cones are located at the corners of 
the zone. In an armchair nanoribbon the condition Eq. (26) 
imposes the quantization of k y . This splits the whole spec- 
trum into a finite number of ID branches. A branch is repre- 
sented by a horizontal line. Gapless branches (solid lines) pass 
through Dirac cones whereas branches with the gap (broken 
lines) do not. A given branch may be represented by more 
than one line on this figure [see discussion before Eq. (32)]. 




FIG. 3: A segment of an armchair nanoribbon of width W. 
The auxiliary sites at the edges (where the wave function must 
vanish) are shown by the hatched circles. The nanoribbon 
unit cell is enclosed inside the dotted line. The deformed 
bonds at the edges are shown bold. 
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FIG. 4: Qualitative behaviour of the conductance of the long 
nanoribbon as a function of disorder (np is the concentra- 
tion of the disordering radical). The dashed line marks the 
crossover from the gap-dominated to the disorder-dominated 
regime. Note the counter- intuitive trend left of the dashed 
line: the conductance is increasing with increasing disorder! 



